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微型 摆 云 发 动机 上 同 隐 内 部 流动 数值 分 析 
孙 萌 ”张震 宇 : 孔 文 俊 
(中国 科 学 院 工程 热 物理 研究 所 ， 北 京 100190; “中 国 科 学 院 大 学 ， 北 京 100039) 

摘要 本 文采 用 离散 速度 方向 模型 对 微型 摆 式 发 动机 摆 臂 与 腔 体 间 隙 内 的 微 尺度 气体 流动 开展 了 数值 研究 ,分 析 了 摆 臂 运动 对 
隙 内 压 差 驱 动 下 气体 流动 阻力 和 流量 的 影响 规律 , 并 探究 了 气体 稀薄 性 的 影响 。 研 究 发 现 , 当 摆 臂 运 动 方向 与 压 降 方向 相同 时 ， 
稀薄 气体 效应 导致 的 速度 滑 移 相 对 于 气体 流动 阻力 对 流量 的 变化 起 主导 作用 ; 当 摆 臂 运动 方向 与 压 降 方 向 相反 时 ， 由 于 微 尺度 
下 壁面 作用 的 增强 ， 在 边界 附近 出 现 气 体 回流 ， 且 在 气流 进出 口 处 出 现 涡 结构 。 
关键 词 间 辽 流动， 离散 速度 方向 模型 ， 稀 薄 效 应 
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On Clearance Gap Flow in a Micro Internal Combustion Swing Engine 
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AbstractNumerical simulations on clearance gap flows between the swing arm and the engine body were completed 
by using the Discrete Velocity Direction(DVD) model in a Micro Internal Combustion Swing Engine. The 
influences of the movement of the swing arm on the gas flow resistance and fluxes driven by the differential pressure 
within the gap wereinvestigated, and the effects of rarefication were discussed. The results show that when the swing 
arm moves in the same direction with the drop of pressure, the velocity slip caused by rarefication plays a leading 
role in the variation of gas flow resistance and fluxes.When the swing arm moves in the opposite direction with the 
drop of pressure,the recirculation flow and vortex occurred near the boundary and entrance respectivelybecause of 
the enhancement of the wall effects. 
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0 引言 
随 着 科技 的 进步 ， 电 子 和 机 械 产品 趋 于 微型 化 ， 和 摩擦 阻力 Bl， 王 文 等 假设 间隙 中 气流 为 伯 努 利 流 ， 
对 能 量 供给 单元 提出 更 高 要 求 。 微 型 摆 式 发 动机 由 ”采用 管道 流动 的 计算 方法 估算 阻力 损失 外; 郭志 平等 


于 结构 简单 、 效 率 高 成 为 最 具 潜 力 的 一 种 微型 能 源 对 采用 FLUENT 软 件 模拟 间隙 内 的 气体 流动 种。 分 析 
动力 系统 。 但 系统 尺度 减 小 导致 动静 部 件 之 间 的 摩 此 前 的 研究 可 以 发 现 ， 对 于 微型 摆 式 发 动机 微 间 院 

擦 和 泄漏 损失 增加 。 当 发 动机 正常 运行 时 ， 间 辽 尺 内 气体 流动 阻力 和 泄漏 的 研究 大 多 采用 求解 
度 控制 在 微米 量 级 ， 此 时 间隙 内 气体 Knudsen 数 在 Navier-Stokes 方 程 方法 进行 粗略 估算 , 很 少 考虑 微 尺 
0.01~1 之 间 ， 处 于 滑 移 流 甚至 过 渡 流 领域 , 稀薄 气体 。 度 气 体 流 动 的 稀薄 效应 ， 而 速度 滑 移 和 压缩 性 增强 


效应 显著 。 即 使 采用 二 阶 或 更 高 阶 滑 移 边界 条 件 ， 等 微观 流动 特性 会 对 发 动机 的 气动 性 能 产生 显著 影 

适用 于 连续 介质 模型 的 Navier-Stokes 方 程 也 仅 能 扩 啊 。 

展 到 Kn<0.5 的 流动 领域 ， 且 误差 较 大 [1。 本 文采 用 精度 更 高 的 动力 学 方法 一 离散 速度 方 
已 经 有 学 者 对 微型 摆 式 发 动机 微 间 隐 内 摩擦 损 向 模型 "研究 微型 摆 式 发 动机 间隙 内 的 气动 特性 ， 


失 和 泄漏 开展 了 研究 。Gu 假 设 间 际 内 气体 流动 处 于 考虑 微 尺度 气体 流动 的 稀 注 效应 ， 探 究 了 摆 辟 运动 
连续 流 领 域 , 用 Navier-Stokes 方 程 估算 了 间 际 泄漏 量 。 及 气体 稀 注 性 对 微 间 隙 内 气体 流动 阻力 和 流量 的 影 
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响 规 律 ， 为 改进 微型 摆 式 发 动机 结构 ， 提 高 气动 效 
率 提供 参考 。 
1 离散 速度 方向 模型 验证 
1.1 离散 速度 方向 模型 简介 
离散 速度 方向 (DVD) 模 型 是 一 种 简化 Boltzmann 
方程 的 动力 学 方法 , 它 通 过 降低 Boltzmann 方程 的 维 
数 来 减 小 方程 数值 求解 的 计算 量 。DVD 模型 在 保持 
分 子 速率 连续 的 前 提 下 ， 将 分 子 运动 方向 离散 ， 用 
离散 方向 上 的 分 子 速 率 代 检 Boltzmann 方程 中 连续 
的 速度 分 布 空间 , 将 六 维 的 Boltzmann 方程 降低 到 三 
维 ， 从 而 减 小 Boltzmann 方程 数值 求解 的 计算 量 。 
1.2 DVD 模型 验证 
摆 式 发 动机 摆 辟 与 生体 间 的 气体 流动 可 近似 认 
为 是 Couette 流动 和 Poiseuille 流动 的 结合 ， 现 分 别 
在 两 种 流动 中 验证 模型 精度 ， 并 将 数值 计算 结果 与 
线性 化 Boltzmann 方程 (LBE) 方 法 对 比 。 计 算 中 分 子 
离散 速率 个 数 均匀 选取 8 个 ， 离 散 速度 方向 选取 
4x8 个 。 用 各 速率 区 间 的 分 子 数 n, 作为 控制 方程 变 
量 可 得 DVD 模型 的 控制 方程 为 
On; 了 On 2 1 
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计算 得 到 的 无 量 纲 剪 切 力 与 LBE 方法 计算 结果 
由 对 比如 图 1 所 示 。 
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图 1 无 量 纲 剪 切 力 随 Knudsen 数 变化 
Fig.1Dimensionless shear stress Variation with Knudsen number 
1.2.2Poiseuille 流动 

计算 得 到 的 无 量 纲 流量 与 LBE 方法 计算 结果 中 
对 比如 图 2 所 示 。 
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图 2 无 量 纲 流 量 随 Knudsen 数 变化 
Fig.2Dimensionless mass flux variation with Knudsen number 
由 图 1、2 可 知 ， 从 滑 移 流 到 过 渡 流 领域 DVD 
模型 计算 得 到 的 无 量 纲 剪 切 力 和 流量 与 LBE 方法 最 
大 误差 均 不 超过 3.2%， 计 算 精 度 明 显 高 于 滑 移 边界 
条 件 的 Navier-Stokes 方程 。 

2 微型 摆 式 发 动机 间隙 内 的 气体 流动 
在 微型 摆 式 发 动机 摆 臂 与 仙人 体 之 间 的 间隙 内 ， 
由 于 摆 臂 运动 和 压 差 同时 存在 ， 其 中 的 气体 流动 属 
于 Couette-Poiseuille 流 。 腔 体内 高 温 燃 气 进 入 间隙 内 ， 
若 间隙 尺寸 为 SU ，Knudsen 数 即 可 达 0.1 左右 ， 若 
间隙 尺寸 进一步 减 小 ， 则 进入 过 渡 流 领域 。 现 分 别 
研究 摆 臂 运动 和 和 气体 稀薄 性 对 流动 阻力 和 流量 的 影 
响 规律 。 以 氮气 为 例 进行 计算 ， 初 始 压力 为 标准 大 


气压 101325Pa ， 温度 为 273.15K ， 压 差 


P3027x107 Pa/m. 
dx 


2.1 摆 尽 运动 对 气体 流动 阻力 和 流量 的 影响 
计算 中 改变 摆 辟 运动 方向 和 速度 ， 得 到 间 际 内 
气体 质量 流量 和 流动 阻力 的 变化 曲线 如 图 3、4 所 示 。 
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3 间隙 内 气体 质量 流量 随 摆 臂 运动 速度 的 变化 


Fig.3 Mass flux variation with the Swing arm movement 
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W/m/s 
4 间隙 内 气体 流动 阻力 随 摆 臂 运动 速度 的 变化 

Fig.4Flow resistance Variation with the Swing arm movement 
由 图 3、4 可 知 ， 在 同一 Knudsen 数 情况 下 ， 无 
论 摆 臂 正 向 或 反 向 运动 ， 间 隙 内 气体 流动 阻力 均 随 
摆 臂 运动 速度 增 大 而 增 大 ， 可 见 压 差 对 流动 阻力 的 
影响 较 小 ; 在 摆 臂 运动 速度 相同 情况 下 ， 随 Knudsen 
数 增 大 间 隐 内 气体 流动 速度 梯度 增 大 导致 流动 阻力 
增 大 。 当 摆 臂 正 向 运动 时 ， 和 气体 流量 随 摆 臂 运 动 速 
度 增 大 而 增 大 ， 当 摆 臂 反 向 运动 时 ， 和 气体 流量 随 摆 
臂 运动 速度 增 大 先 正 向 减 小 后 反 向 增 大 。 和 气体 流量 
随 摆 臂 运动 变化 的 原因 可 由 间隙 内 气体 流 线 图 看 出 。 


图 5Kn=0.1128 摆 臂 反 向 运动 速度 0.1myAs 的 流 线 图 
Fig.S Streamline with v=0.1m/s, Kn=0.1128 


pa 


6Kn=0.1128 摆 臂 反 向 运动 速度 0.15m/s 的 流 线 图 
Fig.6 Streamline with v=0.15m/s, Kn=0.1128 
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图 7Kn=0.3385 摆 臂 反 向 运动 速度 0.1m/s 的 流 线 图 
Fig.7 Streamline with v=0.1m/s, Kn=0.3385 
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图 8 Kn=0.3385 摆 臂 反 向 运动 速度 0.15SmyAs 的 流 线 
Fig.8 Streamline with v=0.15nys, Kn=0.3385 
昌 图 5-8 可 知 , 当 摆 辟 反 向 运动 时 在 近 壁 区 气体 
出 现 回流 。 当 Knudsen 数 相同 时 ， 随 摆 臂 运动 速度 
增 大 回流 区 逐渐 扩大 ; 当 摆 臂 运动 速度 相同 时 ， 随 
Knudsen 数 增 大 回流 区 逐渐 扩大 。Kn=0.1128 时 ， 回 
流 区 扩大 至 主流 区 域 ， 且 在 气流 的 进出 口 处 出 现 明 
显 的 涡 结 构 ， 此 时 壁面 运动 与 压 差 作用 对 气体 流动 
的 驱动 作用 几乎 相同 ; Kn=0.3385 时 ,气体 全 部 反问 
流动 ， 此 时 壁面 作用 的 影响 已 扩大 至 全 部 流 场 区 域 ， 
间 隐 泄漏 量 主要 由 摆 臂 运动 速度 决定 。 由 此 可 知 ， 
对 于 微型 摆 式 发 动机 微 间 隐 内 的 气体 流动 ， 特 征 尺 
度 决 定 了 控制 间 际 泄漏 流动 方式 的 主要 因素 。 
2.2 稀薄 性 对 气体 无 量 纲 流动 阻力 和 流量 的 影响 
计算 中 改变 气体 流动 的 Knudsen 数 ， 和 气体 无 量 
纲 质量 流量 和 流动 阻力 的 变化 曲线 如 图 9、10 所 示 。 
由 图 9 可 知 ， 当 摆 臂 正 向 运动 时 ， 随 气体 稀薄 
性 增强 区 域 中 心 处 气体 速度 迅速 减 小 导致 间隙 内 气 
体 的 无 量 纲 质量 流量 减 小 。 由 图 10 可 知 ， 随 气体 稀 
薄 性 增强 ， 无 量 纲 流动 阻力 增 大 ， 摆 臂 运动 速度 和 
方向 对 流动 阻力 影响 较 小 。 
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Fig.9 Dimensionless mass flux variation with Knudsen number 
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图 10 无 量 纲 流 动 阻力 随 Knudsen 数 的 变化 


Fig.10 Dimensionless shear stress variation with Knudsen 


number 

3 结论 

本 文采 用 离散 速度 方向 模型 研究 微型 摆 式 发 动 
机 摆 辟 与 生体 之 间 微 间 隐 内 的 微 尺度 气体 流动 ， 探 
究 了 气体 稀薄 性 和 摆 臂 运动 对 间隙 内 压 差 驱动 下 的 
气体 流动 阻力 和 流量 的 影响 。 主 要 得 出 如 下 几 点 结 
论 : (]) 摆 臂 运 动 方向 与 压 降 方向 相同 时 ， 随 摆 臂 运 
动 速度 增 大 ， 气 体质 量 流量 增 大 ; (2) 摆 臂 运动 方向 
与 压 降 反 向 相反 时 ， 随 摆 臂 运动 速度 增 大 ， 和 气体 质 
量 流量 先 减 小 后 增 大 ， 近 壁 区 及 进出 口 处 出 现 回流 
区 和 涡 结 构 ; (3) 无 论 摆 臂 运动 方向 与 压 降 方向 是 否 
相同 ， 气 体 流动 阻力 均 随 摆 辟 运动 速度 增 大 而 增 大 。 
对 微型 摆 式 发 动机 间隙 内 的 稀薄 气体 流动 而 言 ， 当 
摆 臂 运动 方向 与 压 降 方向 相同 时 ， 应 尽量 降低 摆 臂 
运动 速度 ， 以 减 小 间 隐 气体 泄漏 量 ， 而 当 摆 臂 运 动 
方向 与 压 降 方向 相反 时 ， 应 综合 考虑 间隙 内 气体 流 
动 阻力 和 泄漏 量 对 整 机 效率 的 影响 ， 设 计 摆 臂 的 最 
佳 运动 速度 以 提高 系统 的 运行 效率 。 
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